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Vector models of gravitational Lorentz symmetry breaking 
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Spontaneous Lorentz symmetry breaking can occur when the dynamics of a tensor field cause it to 
take on a non-zero expectation value in vacuo, thereby providing one or more "preferred directions" 
in spacetime. Couplings between such fields and spacetime curvature will then affect the dynamics 
of the metric, leading to interesting gravitational effects. Bailey & Kostelecky [f ] developed a post- 
Newtonian formalism that, under certain conditions concerning the field's couplings and stress- 
energy, allows for the analysis of gravitational effects in the presence of Lorentz symmetry breaking. 
We perform a systematic survey of vector models of spontaneous Lorentz symmetry breaking. We 
find that a two-parameter class of vector models, those with kinetic terms we call "pseudo-Maxwell," 
can be successfully analyzed under the Bailey-Kostelecky formalism, and that one of these two 
"dimensions" in parameter space has not yet been explored as a possible mechanism of spontaneous 
Lorentz symmetry breaking. 

PACS numbers: 03.50.-z, 04.25.Nx, 04.50.Kd 



I. INTRODUCTION 

It is widely believed that classical general relativity, 
as formulated by Einstein, is a particular limit of some 
underlying theory of quantum gravity. However, at en- 
ergy scales that are now accessible, it is expected (from 
our knowledge of effective field theory) that any funda- 
mentally non-classical effects would be suppressed by at 
least a factor of the ratio of our experimental energy 
scale to the Planck scale; even for today's most pow- 
erful particle colliders, this ratio still gives a suppression 
factor of 10~ 16 . With no foreseeable way to bridge this 
sixteen-order-of-magnitude gap in energy, we are forced 
to aim for sensitivity rather than power when searching 
for quantum-gravitational effects. 

One particularly interesting avenue for this search is 
the possibility of quantum-suppressed Lorentz violation. 
In such a scenario, the underlying theory would include 
a tensor field (or fields) which spontaneously takes on 
a non-zero expectation value. Such a field would, in 
essence, provide a "preferred" direction or directions in 
spacetime. 1 The background value of this field could then 
couple weakly to conventional matter fields [2]; thus, the 
effects of such a tensor field could in principle be seen 
via careful observation of the behaviour of conventional 
particles and fields. 

A particularly interesting venue in which to search for 
possible violations of Lorentz invariance is the gravita- 



tional sector. Interactions between a dynamical metric 
and a tensor field with a non-zero expectation value have 
been postulated as a possible method of modifying cos- 
mology [3, 4, 5, 6, 7], as a mechanism for modifying New- 
tonian gravity to solve the dark-matter problem [8, 9], or 
simply in their own right as modifications of conventional 
gravity [10, 11, 12]. Such modifications of gravity will, 
in general, cause modifications to the weak-field limit of 
gravity. The linearized effects of a direct coupling be- 
tween Lorentz-violating fields and the Riemann tensor 
were analyzed in some detail by Bailey and Kostelecky 
[1]. By making certain assumptions about the properties 
of the equations of motion, they were able to obtain an 
effective linearized gravitational equation of the form 



SG ab + % ab cdef 6R 



cdef 



8irG6T ab 
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where 8G ab and 6~R a bcd are (respectively) the Einstein 
and Riemann tensors linearized about a flat background, 
5T ab is the stress-energy of conventional matter, and 
^ab c is a "small" tensor (in a sense we will make ex- 
plicit below) depending in a particular way on the back- 
ground values of the Lorentz-violating tensors. Using 
this effective equation, they then performed a thorough 
post-Newtonian analysis of such theories, examining the 
effects of Lorentz-violating fields on phenomena includ- 
ing satellite orbits, interferometric gravimetry, torsion- 
balance experiments, and frame-dragging. 2 

While this formalism is highly valuable for the anal- 
ysis of the interface between gravity and Lorentz viola- 
tion, its range of applicability is not immediately clear. 



'Electronic address: mdseifer@indiana.edu 

1 Such a field is often said to be "Lorentz-violating" . This descrip- 
tion plays somewhat fast and loose with usual notions from the 
rest of physics; the field does, after all, transform as a proper 
tensor field under local Lorentz transformations. A more ac- 
curate way to describe such a field would be to say that it 
"spontaneously breaks Lorentz symmetry", but such phrasing 
is rather more awkward. In the interests of readability and con- 
sistency with other papers in the literature, wc will use "Lorentz- 
violating" in this sense as well. 



2 It is important to note that although the Bailcy-Kostelccky for- 
malism can be applied to the analysis of post-Newtonian gravity, 
the theories to which this formalism can be applied are in gen- 
eral not the same as those to which Will's familiar Parametrized 
Post-Newtonian (PPN) formalism [13, 14] can be applied. The 
connections and distinctions between these two formalisms are 
explored in Section III C of Bailey and Kostelecky's original pa- 
per [1], 
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To obtain the effective gravitational equation (1), it was 
necessary for Bailey and Kostelecky to place certain con- 
ditions on the equations of motion, rather than on the 
action from which they were derived. As action princi- 
ples tend to be conceptually simpler than the equations 
of motion derived from them, it would be quite helpful 
to know whether a given action which includes sponta- 
neous Lorentz symmetry breaking is analyzable in the 
Bailey-Kostelecky formalism. Should this be the case, 
the physical predictions of their paper [1] would be di- 
rectly applicable to any such model. 

This question is the focus of the present work. We 
will restrict our attention to the simplest type of tensor 
field which can spontaneously break Lorentz symmetry, 
namely vector fields A a . In Section II, we describe the 
properties of the theories we will be concerned with, and 
we review the conditions required for successful use of the 
Bailey-Kostelecky formalism. Section III is dedicated to 
the application of these conditions to the vector actions 
under consideration; we will see that the class of vector 
theories for which the Bailey-Kostelecky formalism can 
successfully be used is not large, but that there do exist 
previously unconsidered models which can be analyzed 
in this framework. Finally, we discuss these results in 
Section IV. 

We use the sign conventions of Wald [15] throughout, 
and units in which c = 1 . 



II. EQUATIONS OF MOTION AND 
FORMALISM 

A. Actions for Lorentz-breaking vector fields 

Bailey and Kostelecky's analysis of gravitational 
Lorentz violation [1] begins by assuming an action of the 
form 




Ceh here is the usual Einstein- rlilbert action, 

C EH = R - 2A. (3) 

We will assume throughout that A = 0. The second term, 
Clv, contains the non-trivial couplings of the Lorentz- 
violating fields to the metric: 

C LV = -uR + s ab {R T ) ab + t ahcd C abcd . (4) 

Here, R is the Ricci scalar, (R T ) ab = Rat — \g a bR is 
the trace- free Ricci tensor, and C a bcd is the Weyl ten- 
sor. The tensors u, s ab , and t abcd may be fundamental 
fields or (as will be the case in our analysis) composites 
of other fields present in the theory. The final term, C, 
contains the terms determining the dynamics of the fun- 
damental Lorentz-violating fields, as well as the action 
for conventional matter. 



In the case of a single vector field being responsible 
for Lorentz symmetry breaking, we can be more specific 
in the form of the Lagrangian. Denoting the Lorentz- 
breaking vector field by A a , the most general Lorentz- 
violation coupling terms will be of the form 

Clv = i(-fu{A 2 )R + f s (A 2 )A a A»R ab ), (5) 

where A 2 = A a A ai /« and f s are arbitrary functions of 
A 2 , and £ is a coupling constant. 3 (By the symmetries 
of the Wcyl tensor, any term analogous to t abcd C a bcd and 
constructed out of A a and the metric must vanish.) This 
term is best thought of as a "weak" coupling term be- 
tween the vector field and the curvature; the "weakness" 
of this coupling will be of importance in the next subsec- 
tion. 

The dynamics for A a , meanwhile, will be determined 
by C . We can write the Lagrangian for an arbitrary 
second-differential-order vector theory as 

£' = K\ c d V a A b V c A d - V(A 2 ) + 2 K £ mat , (6) 

where £ ma t is the Lagrangian for "conventional" matter; 
k = 8irG; V(A 2 ) is the potential for the vector field, con- 
structed to have a minimum at a non-zero value of A a ; 
and K a b c d is a tensor constructed out of A a and the met- 
ric. This tensor can be taken to be symmetric under the 
simultaneous exchange of a <-> c and b <-> d. The con- 
ventional matter action £ ma t can, in principle, contain 
direct couplings to A a . (We will introduce an explicit 
parametrization for K a b c d in the next subsection.) 

We can easily obtain the Euler-Lagrange equations as- 
sociated with this action by varying the action with re- 
spect to g ab and A a ; there result the equations 

{£ g ) ab = -G ab + U ab + ZB ab + (T A ) ab + «(T mat ) afc = 

(7) 

and 

(£ A ) a = 2H(-fLA a R + f' s A b A c R bc A a + f s A b R ba ) 
+ M b c d ea V b A c V d A e - 2V b (K b a c d V c A d ) 

-2V'A a + 6 -^=0, (8) 

where 

A ab = f u G ab + f' u A a A b R 

+ ±f s g ab A c A d R cd + f s A a A b A c A d R cd , (9) 

- )p{J s A a A b ) + vM a {f s A b ^A c ), (10) 



3 Note that /„ is associated with the Ricci tensor in our 
parametrization, while in Bailey &; Kostelecky's original paper 
the tensor s ab is associated with the trace-free Ricci tensor. 
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{T A ) ah = M c d e f ab V c A d V e Af 

+ V e ((A" c d (a|e| A b) - K c d e{a A b) - K^ c d A e )V c A d \ 



- g ^v - A a A b V\ (11) 



(Tniat ) 



M c 



ab _ 1 1 gCy/fgAnat) 
2 ^ffaft 



ab jv-c e 

2-9 



and 



b d 



5K 



b d 



SA a 



(12) 



(13) 



(14) 



(The arguments of the functions /„, f s , and V will be 
regularly omitted for brevity hereafter.) 



B. Bailey-Kostelecky Formalism 

The basic tack taken by Bailey and Kostelecky in their 
original paper [1] was to start from an action of the form 
(2), with its associated equations of motion; to construct 
the linearized equations of motion about a particular type 
of background; and to then impose certain conditions on 
the background and the equations of motion such that the 
linearized equations could be reduced to a particularly 
simple form: 

SG ab = K,(6T mat )ab + uSGab + riabS cd 5R c d 
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- 2s c {a SR b]c + -s ab 5R + s cd 6R acdbl (15) 

where u and s ab are the background values of the fields 
u and s ab . We now review and discuss these conditions 
as they pertain to the vector theories we are considering. 

1. The background values of the Lorentz-violating 
fields are constant with respect to a background flat 
spacetime. In other words, if e is our linearization 
parameter, we are looking for a family of solutions 
such that 

g a b = Vab + eh ab A a =A a + eA a (16) 

with A a 7^ 0, and, in addition, that 

V a A b ~ 0( e ). (17) 

We will see below that these requirements constrain 
the background values of V, as well as greatly sim- 
plifying the equations of motion (7) and (8). 

2. The dominant Lorentz-violating effects are linear 
in the vacuum values u, s ab , and t abcd . This can 
be enforced in our case by working only to linear 



order in the coupling constant £, discarding terms 
of 0(£ 2 ) or higher. Turning this condition around, 
we will also require that in the limit of vanishing £, 
the metric will obey the Einstein equations; this en- 
sures that our "Lorentz-violating" perturbed met- 
ric will only differ slightly from the usual perturbed 
metric derived from the conventional Einstein equa- 
tions. 

3. The fluctuations u, s ab , and t abcd of the Lorentz- 
violating fields do not couple to the "conventional 
matter" sources. This can be ensured by demand- 
ing that 



5£ r 



5A a 



0. 



(18) 



thereby eliminating the last term from equation (8) 
above. In essence, this requirement ensures that it 
is only the metric that is directly affected by the 
dynamical Lorentz breaking. "Conventional" test 
particles will still move on geodesies with respect to 
the now-distorted metric, and these distorted paths 
can in principle allow us to indirectly observe the 
effects of Lorentz violation on gravity. In the re- 
mainder of this paper, we will be studying "vacuum 
solutions" , with all conventional matter sources set 
to zero. 

The independently conserved piece of the Lorentz- 
violating stress-energy (T A ) ab vanishes. More 
specifically, if we take the divergence of the Einstein 
equation (7), we find that the divergence of (Ta)" 6 
must equal the divergence of £ (A ab + B ab ) . This re- 
lation then allows us to "reverse-engineer" the form 
of (TA) ab , up to a piece S ab whose divergence van- 
ishes. This condition is then the statement that 
S ab itself vanishes. 4 

When the Einstein equation (7) is linearized, any 
second derivatives of A a can be eliminated from B ab 
and (TA) ab in favour of second derivatives of the 
metric. In practise, this elimination can only occur 
via the linearized vector equation of motion. This 
condition will be our primary focus in Section III. 

As a consequence of the first condition above, the back- 
ground (zero-order) equations of motion reduce simply to 



5. 



1 



and 



rj ab V(A 2 ) + A a A b V'(A 2 ) = 



V'{A 2 )A a = 0, 



(19) 



(20) 



4 Note that this is not strictly speaking necessary for the analysis 
performed by Bailey and Kostelecky to still be valid, as noted 
in the original paper; in fact, it does not hold for the bumblebee 
model [1]. 



4 



which together imply (as would be expected) that of motion then becomes 
V{A 2 ) = V'(A 2 ) = 0. The linearized Einstein equation 



SG ab + £ f fJG ab - f' u A a A b 5R + -f sV ab A c A d SR cd + f' s A a A b A c A d SR cd j 

+ (iQR abc d e + QK abc d e )S(\7 e \7 c A d ) - V"A a A b (2A c A c + h cd A c A d ), (21) 



S(£ g r b = 



where 

Q R abc d e ee 2/;(77 a y e - ^rf^Ad + f s (-rf b A^ d 



and 

Q K abc d e = A {a K b ^ ec d - A^ a K^ c d - A e K {ab ^ c d . (23) 

The linearized vector equation of motion, meanwhile, be- 
comes 

\s(£ A ) a = £,(-f u A a 5R + f s A b A c SR bc A a + f s A b SR ba ) 
- K» a c d 5(X7 b X7 c A d ) - V"A a {2A b A b + h bc A b A c ). (24) 

In equations (21)-(24), the arguments of the functions /„, 
f s , and V, as well as the tensor K a b c d , are understood 
to be evaluated at their background values A a — ► A a and 
fjab — ► Vab', indices are raised and lowered by the flat-space 
metric r] ab . The quantity d(\7 a \7 b A c ) is given in terms of 
flat-space derivatives and the metric perturbation h ab by 

s(w a w b A c ) = d a (d b A c + (d {b h d) c - ^d c h bd ^j A d ^j . 

(25) 

Note that by Condition 1 above, this is an 0(e) quantity. 5 
The quantities SR a b, 8G a b, and 5R = r) ab SR ab , finally, 
are the linearized Ricci tensor, Einstein tensor, and Ricci 
scalar associated with the metric perturbation h ab . 

It will be to our advantage to introduce a concrete 
paramctrization for the tensor K a b c d . Any tensor with 
the proper index structure constructed out of A a and the 



It is also important to note that the flat-space derivative oper- 
ator d a and the covariant derivative operator V a differ only at 
order e. In particular, this means that the covariant derivative 
of an 0(e) quantity (such as V a A b ) differs from its flat-space 
coordinate derivative by 0(e 2 ), which for the purposes of this 
paper is negligible. 



- g ce A (a 8 b) d + 77 c(a A b ^ e d + n c{ - a 5^ d A e ^ 

+ f' s {-T 1 ab A c A e - A a A b T 1 ce + 2T 1 c{a A b ^A e )A d (22) 
I 

metric will be of the form 

K\ c d = £ 1 (A 2 )g ac g bd + € 2 (A 2 )S a b S c d + € 3 (A 2 )8 a d S c b 
+ € 4 (A 2 )A a A c g bd + ^€ 5 (A 2 ) (A a A d 5 c b + A c A b S a d ) 

+ € 6 (A 2 )A b A d g ac + ^£ 7 (A 2 ) (A a A b S c d + A c A d S a b ) 

+ <£ 8 (A 2 )A a A c A b A d , (26) 

(This particular parametrization is due to Zlosnik et al. 
[9].) However, due to the geometric identity 

V Q [f(A 2 )(A a \7 b A b - A b X7 b A a )] 

= f(A 2 ) ((\7 a A a ) 2 - \7 b A a \7 a A b - R ab A a A b ) 

+ 2f'(A 2 ) (A a A c \7 a A c X7 b A b - A a A c \7 a A b V b A c ) , 

(27) 

we can always eliminate one of £2, £3, £5, or £7 via 
an integration by parts (thereby changing f s as well.) 
Hereafter we will take £2 to vanish. The arguments of 
£i(A 2 ) will also generally be omitted for brevity. 

C. "Pseudo-Maxwell" kinetic terms 

Finally, we note two important properties of the vector 
equation of motion (8) for certain choices of K a b c d . Con- 
sider a kinetic term for which K l - ab * >c d = 0. This places 
restrictions on the £i functions: 

£1 + £3 = (28a) 
£4 = -^£5 = £ 6 (28b) 
£7 = £ 8 = (28c) 

Alternately, this condition implies a kinetic term that can 
be written in the form 

K a b c d V a A b V c A d 
= ±(U l9 ac + U 2 A a A c ){U l9 bd + H 2 A b A d )F ab F cd (29) 
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where F ab = 2V [a A b] , d = ±H\, and U = ±H{H 2 - 
(The signs here are determined by the overall sign of 
£i.) As this kinetic term is simply the familiar Maxwell 
field strength tensor contracted twice with a "generalized 
metric" Hig ab + Ti.2A a A b , we will call such kinetic terms 
(and theories containing them) "pseudo-Maxwell." 

Taking the divergence of the vector equation of motion 
(8) for a general K a b c d and linearizing about our chosen 
background, we find that 

a-fLA a VJR + f' s A»A c A a V a 5R bc + f s A b W a SR ba ) 

- K bac d 5(V a V b V c A d ) - 2V"(A 2 )A a A b 5(V a A b ) = 

(30) 

For an arbitrary vector field A a and an arbitrary metric, 
we know that 

V a V b V c A d = V (o V 6) (V c A d ) 

+ \ {Rabc e V e A d - R abe d W c A e ) . (31) 

It can be then be seen that in the case K^ ba ^ c d = 0, to 
linear order in e the divergence of the vector equation of 
motion is simply 

a-fLA a VJR + f' s A»A c A a V a 5R bc + f s A b W a SR ba ) 

= 2V"(A 2 )A a A b S(V a A b ) (32) 

(note that the quantity in brackets in equation (31) is 
0(e 2 ).) Using the linearized contracted Bianchi identity 
V a 8R ab = 7}V b SR, this last equation is equivalent to 

A a V a (t + \f)j 5R + U' s A b A c 5R bc 

- V"(A 2 )6(A 2 )^J = 0. (33) 

where S(A 2 ) = S(A a A a ) = 2A a A a + h ab A a A b . 

This implies that in the case where K^ c d = 0, if the 
linearised quantity in brackets above vanishes on some 
hypersurface to which A a is non-tangent, this quantity 
will vanish throughout spacetime. (Recall that A a is a 
constant vector field in Minkowski space.) Thus, via an 
appropriate choice of boundary conditions, we can im- 
pose 

& = e (-/: + \f)j sr + u' s A b A c 8R bc 

- V"(A 2 )S(A 2 ) = (34) 

everywhere. 6 This equation can be interpreted as telling 
us how much the vector field moves "up" its potential 



We have abused notation somewhat here, inasmuch as the quan- 
tity bT denned by (34) is not obtained as the linearized varia- 
tion of some quantity T . Nevertheless, we will continue to use 
&T throughout as a reminder that equations involving it are not 
exact, but only hold to linear order. 



(recall that the value of the potential V only depends on 
A 2 ), and so we will call the equation (34) the "massive- 
mode" condition. When combined with the linearized 
vector equations of motion (24), this yields 

\8(Pa)« = UsA b 8G ab - K b a c d 5{V b V c A d ) = 0. (35) 

This massive-mode condition can then be used to im- 
pose further conditions on A a and its derivatives. It can 
be shown (sec Appendix A) that by taking the appro- 
priate combinations of the derivatives of the equation of 
motion, we arrive at the equation 

D a b [A c 6(W [b A c] )} = U s A b A c d [a 5G b]c (36) 

where D a b is the flat-space linear second-order differen- 
tial operator 

D a b = c^a 

+ U {S a b A c A d d c d d + A 2 d a d b - A a A c d c d b ) . (37) 

Thus, the operator O a b applied to the one-form v a = 
A b 6(\7[ a A b i) yields a quantity of order £. The properties 
of Da* (see Appendix A) allow us to conclude that under 
the imposition of appropriate boundary conditions, the 
quantity v a will itself be of order £ as long as 

+ A 2 U) > 0. (38) 

Since wc also have 

A a 5(V b A a ) = \&{V b A 2 ) ~ O(0 (39) 

from the massive-mode condition (34) above, we can con- 
clude that under these assumptions, the quantity 

A a 8{V a A b ) = -2v b + U(V b A 2 ) ~ O(0 (40) 

as well. This condition, along with the massive-mode 
condition (34), will become important in our analysis of 
the effective gravitational equations below. 

III. CONDITIONS ON VECTOR DYNAMICS 
A. The Einstein limit 

1. General case 

Recall the second of Bailey and Kostelecky's conditions 
above: namely, that any Lorentz-violating corrections to 
the linearized Einstein equation are linear in the param- 
eter £. This implies that in the limit £ — > 0, the equa- 
tions of motion (21) and (24) must together imply that 
the conventional linearized Einstein equation is satisfied, 
i.e., that SG ab = 0. In this limit, the equations of motion 
become 

-8G ab +Q K ahc d e 5(V e V c A d )-V"A a A h 5(A 2 ) = 0, (41) 
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with Qk defined as in (23), and 

- K b a c d 5(X7 b X7 c A d ) - V"A a S(A 2 ) = 0. (42) 

We will further allow the functions <Li(A 2 ) to be depen- 
dent on £, defining functions Ci(A 2 ) and T>i{A 2 ) such 
that 

ti = Ci + £Di + 0(£ 2 ). (43) 

For the two equations (41) and (42) to imply the valid- 
ity of the conventional linearized Einstein equation, we 
must be able to eliminate the terms containing second 
derivatives of the vector field from (41) using the vector 
equation of motion (42). Since this must occur for an 
arbitrary perturbation of the vector field, with arbitrary 
derivatives, we conclude that this will only occur if for 
some tensor T ab f , 

Q K abc d e = T ab fK e f c d (44) 

in the limit £ — > 0. If this relation holds, then we can 
combine the linearized Einstein equation and the lin- 
earized vector equation of motion to obtain 

8G ab = ~V"{A 2 )(A a A b + T abc A c )8{A 2 ). (45) 

This further implies that if the conventional Einstein 
equation is to hold in the limit £ — > 0, we must cither 
have A a A b + T abc A c = or 8 (A 2 ) = in this limit. 

What form must this tensor T abc have? For later con- 
venience, we will split it up into pieces of 0(£ ) and 

o(ey. 

T abc = j-abc + ^fabc_ ^ 

Moreover, since we are only concerned with the linearized 
equations, we can take T abc to be composed solely of 
background quantities. Since the only two geometric ob- 
jects "in play" in the background are the vector field 
A a and the flat metric r/ ab , and given the symmetry 
q-abc = rybac m h ercilt m t \ je definition of T abc , we con- 
clude that T abc must be of the form 

T abc = U 1 r 1 ab A c + £/ 2 A ( V )c + U 3 A a A b A c (47) 

and 

f abc = Vm ab A c + V 2 A( a rf> )c + V 3 A a A b A c 7 (48) 

where the coefficients Ui and Vi can in principle be func- 
tions of A 2 . Assuming that 5(A 2 ) ^ 0, the constraint 
that A a A b + T abc A c vanish yields: 

Ux = and U 2 + U 3 A 2 + 1 = 0. (49) 

The question now becomes what form K a b c d can have 
and still satisfy the condition (44). As with our other 
quantities, we will split K a b c d into and O^ 1 ) 

parts: 

K\ c d = (K ) a b c d + iK\ c d + 0{e). (50) 



Note that due to the decomposition (43), (Ko) a b c d or 
K a b c d can be obtained by taking the original definition 
(26) of K a t c d and replacing by Ci or 2?j, respectively. 
Similarly, we will define 

Q K abc d e - (QK ) abc d e + iQ K abc d e + 0(i 2 )- (51) 

In the limit £ = 0, we thus have the condition 

(Q K o) abc d e = T abf (K y f c d (52) 

Both sides of this equation consist of various five-index 
tensors constructed from A a and the metric, with vari- 
ous coefficients given in terms of U 2 and the Ci functions. 
(Their exact forms are given in Appendix B, Equations 
(Bl) and (B2).) Matching these coefficients, we obtain 
a set of eleven equations which the C, functions and U 2 
must satisfy. (We of course want a non-trivial solution for 
the Ci coefficients.) Examination of the resulting equa- 
tions shows that we must have U 2 = —2 and U 3 = AT 2 , 
and that the functions Ci must satisfy 

Ci = -C 3 = -A 2 C 4 = 7;A 2 C 5 and C 7 = (53) 

with Cg and C§ arbitrary. This implies a vector kinetic 
term that can be rewritten in the form 

K\ c d V a A b V c A d 

= Q x {g ac - A- 2 A a A c )(g bd - A~ 2 A b A d )F ab F cd 

+ {G 2 g ab + G 3 A a A b ) V a {A 2 )V b {A 2 ) (54) 

where F ab = 2V[ a Af,j and the coefficients Qi are functions 
of A 2 , related to the Ci functions by C\ = 2C?i , Cq = 
4g 2 -2A~ 2 g 1 , and C 8 = 4Q 3 . 

2. Pseudo-Maxwell dynamics 

In the previous subsection, we assumed that a gen- 
eral form for K a b c d . However, as was noted at the end 
of Section II C, a "pseudo-Maxwell" vector kinetic term, 
satisfying K^ ab ^ c d — 0, will behave somewhat differently. 
The linearized solutions obtained from such an action 
will, with the imposition of appropriate boundary con- 
ditions, also meet additional self-consistency conditions 
due to properties of the linearized equations of motion. 
In particular, in the £ — + limit, the condition (34) be- 
comes 

V"{A 2 )8{A 2 ) = 0. (55) 

This allows us to ignore the constraints (49) on T abc , 
as they were imposed by the requirement that the right- 
hand side of Equation (45) vanish. We therefore only 
have the requirement that the second derivatives of A a 
vanish, as expressed by (44), in order to obtain a valid 
Einstein limit. In this case, the full tensors are given by 
Equations (B3) and (B4) in Appendix B. Once again, we 
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perform the matching of coefficients between these two 
tensors, yielding a set of equations that must be satisfied 
by the Ci and Ui functions. Assuming that C\ ^ —A 2 C4, 
these two tensors will be equal if and only if f/ 2 = —2 and 
Ui = Us = 0. 7 We have thus found two possible vector 
field kinetic terms, given by (29) and (54), for which the 
conventional Einstein limit is recovered in the limit of no 
direct coupling to curvature. 

B. Adding Lorentz violation 

In the above section, we obtained vector actions which 
satisfied Condition 2 above; namely, in the limit of no 
direct coupling to curvature, these actions yielded lin- 
earized equations of motion that implied the conventional 
linearized Einstein equation 8G ab = 0. We now wish to 
"turn on" direct coupling between the curvature and the 
vector field by setting £ ^ and place further constraints 
on the form of these actions. 

Although Condition 2 does not yield any constraints 
on the form of the equations of motion at 0(£), we can 
still constrain the vector action by imposing Condition 
5: we must be able to eliminate the derivatives of A a 
from the metric equation of motion (21) via use of the 
vector equation of motion (24). In particular, the terms 
in (21) which contain derivatives of the vector field can 
be written in the form 

(QK abc d e + tQR ahc d e )5(\/ e \7 c A d ) 

= ((QKv) abc d e + ttQK abc d e + Q R abc d e )) <5(V e V c A rf ) 

(56) 

Using the vector equation of motion (21) and the condi- 
tion (44), we can rewrite this as 

(Q K abc d e + ZQ R abc d e )S(VeVcA d ) 

~ £ (-T abf K e f c d + Q K abc d e + Qn abc d e ) S(\7 e \7 c A d ) 

(57) 

where the "~" symbol here means "up to terms not in- 
volving derivatives of A a ." We can further simplify this 
expression by noting that in an arbitrary spacetime, 

V a V b A c = V {a V b) A c - ^R abd c A d (58) 

or, in our case, 

S(V a V h A c ) = <5(V (a V b) A c ) - ^6R abd c A d (59) 

up to linear order in e. Thus, at 0(£) we only need to 
eliminate the symmetrized second derivatives from the 



7 Note that the case where Ci = —A 2 d is a special case of the 
kinetic term (54) derived in the previous section. 



metric equation of motion (21); the antisymmetrized sec- 
ond derivatives will merely result in contractions of A a 
with the linearized Riemann tensor, which are expected 
if the effective linearized gravitational equation is to be 
of the form (1). This will occur if t abc (the £>(£) contri- 
bution to T abc defined in (46)) satisfies the equation 

Q R ab{c d e) + Q K ab(c d e) = t abf K^ f c ) d + T^ bf k^^ d . 

(60) 

This equation is essentially the O (£) analog of Equation 
(44). 

We can now proceed with the analysis of this equation 
as we did in the £ = limit: we write out the left-hand 
and right-hand sides in terms of various five-index ten- 
sors constructed from rf h and A a , and match coefficients 
to determine the possible forms of the T>i 's and their cor- 
responding T abc tensors. Expressions for the resulting 
tensors are given in Appendix B; the left-hand side of 
(60) is given by equation (B5), while the right-hand side 
is given by (B6). 

1. General case 

In the case where K^ c d £ 0, we found in Section 
III A 1 that the kinetic terms for the vector must be given 
by (54), with U3 = A~ 2 . We now wish to match the 
coefficients in (B5) and (B6) to see what conditions can 
be placed on the T>i coefficients and the functions f s and 
f u . Substituting in the appropriate relations for the C^'s 
and U3, we find that if (B5) and (B6) are to agree, we 
are forced to set 

fL(A 2 ) = (61) 

and 

f s (A 2 ) = 0. (62) 

These conditions can most easily be seen from the coeffi- 
cients of r] e ( a r] b * >c A d and ri ab A^ c S e ^ d , respectively. In other 
words, the vector model whose kinetic term is given by 
(54) cannot be modified with a Lorcntz-violating curva- 
ture coupling of the form (5) and still satisfy the assump- 
tions of the Bailey-Kostelecky formalism. (Note that set- 
ting f u (A 2 ) to a non-zero constant merely changes the ef- 
fective value of G.) Thus, this theory cannot be success- 
fully be analyzed under this formalism unless Lorentz- 
violating effects induced by the coupling term Clv van- 
ish. 



2. Pseudo-Maxwell dynamics 

The obvious next step is to attempt the same coef- 
ficient matching for pseudo-Maxwell vector theories, as 
defined in (29). However, when we naively do so, we find 
that the same logic that forced us to abandon Lorentz vi- 
olation in the vector model (54) again forces the Lorentz- 
violating functions f u and f s to vanish in the case of 
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pseudo-Maxwell kinetic terms. This stands in opposition 
to the fact Bailey and Kostelecky successfully applied 
their formalism to the so-called "bumblebee model" [2] 
in their original paper [1] ; the kinetic term for this model 
is the same as our pseudo-Maxwell kinetic term in the 
special case Ci = constant and £4 = 0. What have we 
failed to take into account? 

The missing pieces are the conditions on the linearized 
derivatives of A a derived in Section II C. Namely we 
found that under the imposition of certain boundary con- 
ditions, we have 

A a S(V a A b ) ~ A a 5(V b A a ) ~ O(0 (63) 

everywhere in the spacetime. The role of these condi- 
tions is easiest to see by returning to Equation (57) and 
examining the C(^ 1 ) derivative terms remaining in the 
equations of motion after eliminating the 0(£°) deriva- 
tive terms. To wit, suppose there exist tensors T ab ^ 
c abc d , c abc d, and d abcd such that we can write 

- T abf K e f c d + Q K abc d e + Q R abc d e 

= f abf (K ) e f c d + c abc d A e + c abe d A c + d abce A d (64) 

The conditions (63) on the derivatives of A a imply that 
to linear order in e, A c S{V e V c A d ) and A d S(V e V ' c A d ) are 
of order £; similarly, to this order in e wc will have 

A e S(V e V c A d ) = A e 5(V c V e A d ) + A e 6R cef d A f 

= A e A f 6R cef d + O(0- (65) 

Thus, if Equation (64) holds, we will have 

£(-% abf k e f c d + Q K ahc d e + QR abc d e )S(\7 e \7 c A d ) 
= e(T a ^(X ) e / c d( 5(V e V c A d ) + c abc d A e AfR cef d ) 

+ 0(e), (66) 

since all the other terms on the right-hand side of (64) are 
of 0(£) when contracted with <5(V eV c A d ). 8 In essence, 
the derivative conditions (63) allow us to "ignore" certain 
of the equations arising from the coefficient-matching im- 
plicit in (60) at a given order in £. 

To perform this decomposition, we first note that by 
taking the equation T^ b ' (K ) e f c d = (Q,Ko) abc d e and re- 
placing the Ci functions with T>i functions, we obtain 

T abf K e f c d = Q K abc d e - (67) 

(To put this another way, the relations (28) hold to all 
orders in £, and so T$ b ^ K e f d = Q,K abc d e to all orders.) 



Note that the decomposition in (64) is ambiguous: it does not 
address what is to be done with terms of the form C abc A d A e , for 
instance. However, it is easily seen from (66) that such terms will 
vanish when contracted with the Ricmann tensor, so it does not 
matter whether we consider them to be part of c a6c £ j or d abce . 



Thus, the first two terms on the left-hand side of (64) 
cancel, and we merely need to examine Qr d e to find 
out the required form of the tensors on the right-hand 
side. The form of Q R abc d e is given by (22); for a f ab f 
given by (48), the quantity T ab f K e f c d is given by 

f abf K e f c d = V 2 C l {A {a 5 b ^ d T 1 ce - A^ a rp c 8 e d ) 
+ ^(Ci + A 2 C 4 )(?7 a y e Ai - v ab S e d A c ) 
+ (V 2 C 4 + Vsid + A 2 C i ))A a A b {r 1 ce A d - A c S e d ) 

+ V 2 C 4 (A^ a S b) d A c A e - A^r 1 b '> c A d A e ). (68) 

Comparing these equations, we can then see that Equa- 
tion (64) is satisfied if f abc has 

V 2 d = -f„ (69) 

with V\ and V3 arbitrary, and 

c abc d = U {-\rf h ^d + V c{a 5 b) d^j ■ (70) 

Note that this latter quantity is independent of the form 
of f ab f. 

Finally, we confirm that the effective gravitational 
equations are of the proper form for these pseudo- 
Maxwell models. Applying the massive-mode condition 
(34) to the linearized Einstein equation (21), we obtain 

6G ab = £ ^f u 5G ab - i/ s (A a A b 5R - rf b A c A d 5R cd )^j 
+ (Q K abc d e + ZQR abc d e )5{V e V c A d ). (71) 

Using the linearized vector equation of motion (35) con- 
tracted with 7^ a6 ^ = —2A ( - a r] b ^, we can eliminate the 
0(£°) derivative terms to obtain 

5G ab = ^fu5G ab + ^f s A a A b 5R - 2f s A^ a 5R b \A c 

+ ^f sV ab A c A d SR cd + Q R abc d e S(\7 e \7 c A d ^j . (72) 

Lastly, the remaining derivatives of A a in the above equa- 
tion can be eliminated using the derivative conditions, as 
noted above in equation (66); this yields 

^- { (,^ + /.(i**M-2*-«M- 

+ if b A c A d 5R cd + A c A d 6R a cd b ^j ^ . (73) 

In our parametrization, the bumblebee model [2] is ob- 
tained by setting / s = 1 and /„ = 0. Plugging in these 
values, this effective equation for SG ab reduces to the 
form of the effective gravitational equation (15) found by 
Bailey and Kostelecky, with an "effective u" of ~jA 2 
and with s ab = A a A b - \if b A 2 . 
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IV. DISCUSSION 

We have systematically examined the dynamics of 
vector-tensor gravity theories with spontaneous Lorentz 
symmetry breaking. The primary constraints on the 
form of these theories were obtained by imposing two 
of Bailey & Kostelecky's conditions: First, we required 
that the equations have the correct weak-field Einstein 
limit SG a b = when the Lorentz-violating terms (5) are 
"turned off" (Condition 2 of the list in Section II B); sec- 
ond, we required that the linearized stress-energy of the 
vector field vanish automatically when the linearized vec- 
tor equations of motion held (Condition 5). The first of 
these requirements led us to the conclusion that the ki- 
netic terms for our vector fields must be of the form (29) 
or (54). The vanishing of the linearized vector stress- 
energy was found to be a somewhat more subtle issue; we 
found that under the imposition of appropriate boundary 
conditions, the so-called pseudo-Maxwell vector models 
(those with kinetic terms of the form (29)) could lead to 
effective gravitational equations expressed solely in terms 
of the metric. 

It is important to reiterate that the imposition of 
boundary conditions is necessary to obtain effective grav- 
itational equations of the form used by Bailey and Kost- 
elecky in their post-Newtonian analysis; as was noted at 
the beginning of Section IIIB2, an arbitrary solution of 
the vector equations of motion will not have the proper 
relations between the derivatives of the vector field to 
cause the linearized vector stress-energy to vanish. In 
a certain sense, this confirms the aptness of the name 
"bumblebee model" . This name was originally inspired 
by the notion that according to received wisdom, bum- 
blebees should not be able to fly; naive calculations by 
engineers and entomologists in the 1930s seemed to show 
that the bumblebee's wings were too small to allow it 
to fly, and only once more subtle aerodynamic effects 
were taken into account was the mystery explained. Sim- 
ilarly, a naive comparison of the bumblebee vector equa- 
tions of motion with its stress-energy causes us to con- 
clude that we cannot introduce Lorentz-violating gravi- 
tational effects into the model; only once more subtle ef- 
fects (namely, proper boundary conditions) are taken into 
account can Lorentz violation in the bumblebee model 
"fly." 

This said, the technique of imposing boundary condi- 
tions to obtain the desired effective gravitational equa- 
tions is not entirely rigourous. In particular, we used 
the somewhat vague statement that "solutions depend 
continuously on initial data" to argue that the quantity 
A a 5(V a Ab) was of order £. While this is true, the notion 
of continuity associated with well-posedness of an initial 
value problem is defined in terms of the norms of the so- 
lutions on certain Sobolev spaces, and is not easy to gain 
a simple intuition about (see Chapter 10 of [15]). The 
notion of "continuous dependence on initial data" (and, 
by Duhamel's principle, on sources) does allow us to say 
that we can always make A a S(V a Ab) as small as we like 



by tuning £ to be "sufficiently small"; however, it is far 
from clear how small is "sufficient." It would be instruc- 
tive to obtain more careful estimates of how critically the 
magnitude of A a 8(V a A\,) depends on £; however, such an 
analysis is well outside the scope of this paper. 

In some sense, the fact that only pseudo-Maxwell ki- 
netic terms are acceptable for Lorentz violation is not 
entirely surprising given the Bailey-Kostelecky formal- 
ism's requirement of cancellations in the equations of mo- 
tion. The quantity V a Ab will, in general, depend both 
on derivatives of the vector field and derivatives of the 
metric (this latter dependence can be thought of as aris- 
ing from the Christoffcl symbols implicit in W a Af,.) A 
vector kinetic term containing an arbitrary contraction 
of W a Ab with itself and other fields will then, in general, 
lead to a "cross term" between derivatives of the vec- 
tor and derivatives of the metric in the kinetic terms of 
the theory [16]. However, the antisymmetrized derivative 
V[ a Ab] is independent of the metric, and so the kinetic 
terms for the metric and the vector will be decoupled 
when we contract V[ Q Ab] with itself. It is therefore not 
surprising that this special property should have some 
bearing on the relation between the vector equations of 
motion and the gravitational equations of motion. 

In the case of €4 = and €1 constant, the pseudo- 
Maxwell theories we have been discussing become a sim- 
ple Maxwell action for the vector field (albeit without 
gauge symmetry, which is broken by the presence of the 
potential.) However, the theories for which (£4 ^ do not 
appear to have been previously considered in the litera- 
ture, at least as far as concerns Lorentz-violating effects. 
In some sense, the presence of a £4 ^ term causes 
Lorentz violation for the Lorentz-violating field itself: at 
the linearized level, small perturbations of the vector field 
"see" the effective metric Hig ab + H 2 A a A b (as defined in 
(29)), rather than the spacetime metric g ab . In particu- 
lar, in the bumblebee model the Nambu-Goldstone modes 
of the Lorentz-violating vector field can be interpreted as 
a Maxwell field in a particular gauge [17]. If we naively 
extended this interpretation to a general pseudo-Maxwell 
theory, one would expect that the "speed of light" would 
be different from the "speed of gravity" , as the two fields 
would propagate on the null cones of two different met- 
rics. Under such an interpretation the "photon" would 
almost certainly propagate anisotropically; it is also pos- 
sible that such an interpretation would predict vacuum 
birefringence. Experimental bounds on such phenomena 
could then place bounds on the relative values of Hi and 
H.2- That said, this intuitional understanding may be 
complicated by the fact that the correspondence in the 
above-mentioned work [17] is in a non-standard gauge. 
It is also known that this correspondence does not carry 
over to theories with more general kinetic terms than the 
bumblebee model [18], though the class of models exam- 
ined in this last work did not include the pseudo-Maxwell 
theories we have found. More work is needed to elucidate 
the correspondence (if any) between Maxwell theory and 
the Nambu-Goldstone modes of these new theories. 
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Finally, it is important to note that our results imply 
that the Bailey-Kostelecky formalism cannot successfully 
analyze theories with non-standard kinetic terms [4, 5, 
8, 11]. This does not imply that post-Newtonian effects 
in such theories cannot be analyzed; in fact, Bailey and 
Kostelecky did precisely this in their original paper [1] 
for a Lagrangian identical to what Carroll et al. later 
called sigma-sether theory [4]. It is further possible that 
such a theory might in fact provide a viable model of 
Lorentz violation, consistent with current experimental 
constraints, even though it does not fit into the Bailey- 
Kostelecky formalism. In the absence of a more general 



J 



formalism for gravitational Lorentz violation, however, 
such theories will have to be analyzed on a case-by-case 
basis. 
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APPENDIX A: DERIVATION AND HYPERBOLICITY OF THE OPERATOR D D 

Consider the following linearized combination of the vector equations of motion: 



\j (A b (V b (£ A ) a - W a (£ A ) b )) = -AfK dec [a 6 b f] 6(V b V c V d A e ) - V"(A 2 )A b A [a W b] S(A 2 ) 



+ £ {-f u A h A {< y b] 5R + f s A b A c A d A [a V b] SR cd - f s A b A c V [a 5R b]c ) 
Writing out the term — A? K dec [ a 5 b f] for a theory in which K^ ab ^ c d = 0, we find 



(Al) 



_ A f K dec [a S b f] = £i(3A[ fc ^ a ?fl c + s^ d a A e ^] bc ) + U{A b A c 5 [d a A^ + A 2 5 b a ^ d A e ^ - A a A b rj c ^ d A^) (A2) 

Since 5(V[hV c ] Vd^4 e ) ~ 0(e 2 ) and V[ a V(,A c i = 0, we can rewrite the first term on the right-hand side of (Al) (to 
linear order) as 



- 2^if de W(V fc V c V d A e ) = -2A^ de VV(V c V fc V^ e ) = 2D a d [A e S(V [d A e] )} 

I 



(A3) 



Further, applying the massive-mode condition 5 J- = 0, 
we can eliminate the term proportional to V"{A 2 ) from 
(Al), yielding 



D a b [A c 6(V [b A c] )} = UsA b A c d [a 8G b]c 



(A4) 



when the linearized vector equation of motion is satisfied. 

Thus, the quantity v a = A b S(V[ a A b -i) will satisfy a 
second-order differential equation (A4) in flat spacetime. 
Moreover, the source for this equation is "small", i.e., of 
order £. We are thus led to the following question: under 
what conditions will the solution for v a itself be of order 
£? More precisely, let us pick some time coordinate t 
on Minkowski space. We know that if we set £ = 0, 
v a = for all t is a valid solution of the Cauchy problem 
for (A4) with the boundary condition v a (to) = and 
dv a /dt\ to = 0. We wish to know whether, as we "tune" 
£ to zero, the solutions of v a go "smoothly" to zero for 
these boundary conditions. 

This is precisely the question of whether the operator 
D a b has a well-posed initial-value formulation. 9 While 



the general problem of whether an arbitrary operator 
possesses an initial-value formulation can be quite subtle, 
for operators in flat spacetime with constant coefficients 
(such as D a b ) the situation is more clear-cut. Suppose 
O a b is a linear m th -order differential operator which op- 
erates on A^-tuplcs of functions in flat spacetime. (Thus, 
an equation of the form O a b v b = is a system of N linear 
?n th -order differential equations.) Associated with any 
such operator we can find an N x N polynomial-valued 
matrix P a b (X,() such that 



= £>„ 



(A5) 



i.e., if we take P b and replace A by d/dt and ( by V, we 
obtain the operator O a b - We will further assume that the 
matrix P a b is constant with respect to space and time. 
It can then be shown [19, 20] that such an operator has 
a well-posed initial value formulation (with respect to 
an initial-data surface t = constant) if and only if there 



9 Note that a "small" variation in the source terms in (A4) can be 



mapped to a "small" 
Duhamcl's principle. 



variation in the boundary conditions via 
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exists a real number c such that the mN roots Xi of the 
equation 



det 



P(i\,i£) 



= 



(A6) 



satisfy 3(Aj) > — c for all real vectors £. Such an operator 
is said to be "hyperbolic in the sense of Garding." 

To apply this result to the case of the operator D a b , let 
us choose a Cartesian coordinate system on flat spacetime 
{t, x, y, z] for which A x = A y = 0. Then the polynomial 
defined by (A6) becomes 



(g: 1 + a 2 £ 4 )(a^c 



(A7) 



This polynomial has roots when X 2 = ( 2 due to its sec- 
ond factor; these will obviously have 9(A;) = for all 
real (. The third factor, meanwhile, is a slightly more 
complicated quadratic polynomial in A; its roots can be 
shown to be real if its discriminant is positive: 

£ = d ((Ci - U (A 4 ) 2 ) Cl + (£i + UA 2 ) C) > 0, 

(A8) 

where Ci = C 2 + C 2 - If the quantity J) is negative for 
some value of C, the imaginary part of these roots will be 
iVSX Moreover, should this quantity D be negative for 

I 



some real vector C, the magnitude of the imaginary part 
of these roots can be made arbitrarily large: if $i(Ai) = 
±V®o for a given <f = Co, then 3(A t ) = ±Mv^ for 
C = M(o- Thus, the operator D a b defined in (37) will be 
hyperbolic in the sense of Garding if and only if 2) is a 
positive definite quadratic form in C, i.e., if 

- U {A 1 ) 2 ) > and + UA 2 ) > 0. (A9) 

We can therefore conclude that in any frame in which 
these inequalities hold, we can then impose boundary 
conditions on some initial-time surface t = to such that 
^4 b ^(V[ a A b ]) ~ C(C) throughout the spacetime. We can 
further ask that such a frame have A t ^ 0; if this is the 
case, then the massive-mode condition (34) can also be 
imposed on the surface t = to, and it will follow (via the 
linearized equations of motion) that the massive-mode 
condition is satisfied everywhere. Such a frame will nec- 
essarily exist if 



Ci(Ci + A 2 U) > 0. 



(A10) 



(If A 2 < 0, the frame in which A z = satisfies our 
requirements: if A 2 > 0. the required frame is one in 
which A* is non-zero but sufficiently small that <t 2 > 
<ti<ti{A 1 ) 2 .) For A 2 ^ 0, this is equivalent to the con- 
dition that the "effective metric" appearing in (29) is of 
signature (— + ++) or (+ — ). 



APPENDIX B: TENSOR COEFFICIENT-MATCHING 



For a general vector theory, we will have 

(Q K o) abc d e = (Ci - C 3 )A^ c S e d + (C 3 - d)I^ 6 V e - (d + C 3 )5^ d ^ c A e 

- t£b) A a A b A c S e d - 2C 4 A^ a S b) d A c A e + fic 8 - C 6 ) A a A b rf e A d 



~ CsA^rj^AdA* - -C 7 A a A b S c d A e - -C 7 r) ab A c A d A e - C s A a A b A c A d A e . (Bl) 

Assuming that K^ ab ' )c d ^ 0, the tensor T abc must have Ui = and U 2 + U3A 2 + 1 = 0; multiplying these two tensors 
together, we find that 

T abf (K ) e f c d = U 2 C 3 A^ c S e d + UtCxA^hT ~ (\c 5 + i~ 2 (l + U 2 )C 3 ) A a A b A c S e d 



hz 7 A a A b 5 c d A e + U 2 C i A { - a 5 b ^ d A c A e - (C e + A- 2 (l + U 2 )C 1 ) A a A b v ce A d + ^U 2 C 5 A^ a ^ c A d A e 

+ hj 2 C 7 A^if^A c A d - (Cs + A- 2 (l + U 2 ) \C 4 + l(C 5 +C 7 )X\ A a A b A c A d A e . (B2) 



For a pseudo-Maxwell vector theory, we can obtain (Q,Ko) abc d e simply by applying the conditions (28) to (Bl); the 
result is 

{Q K oT hc d e = 2C 1 A( a 7 1 b '> c 5 e d - 2C 1 A {a 5 b ^ d r ] ce + 2C 4 A a A b A c S e d - 2C i A ( ~ a 8 b) d A c A e - 2C 4 A a A b T 1 ce A d + 2C i A {a ^ c A d A e . 

(B3) 
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Due to the massive-mode condition, however, the above constraints on the functions Ui are relaxed; we thus must 
allow for arbitrary C/j functions, yielding 

T abf (K ) e f c d = U 2 C 1 (A^ d r ] ce - A^ c S e d ) + U x {& + A 2 C 4 ) v ab ( v ce A d - A c 6 e d ) 

+ (U 2 C 4 + U 3 (d + A 2 C 4 )) A a A b (rf e A d - A c S e d ) + U 2 C 4 {A^ a 8 b \A c A e - A^r) b >A d A e ). (B4) 

At 0(£), we can attempt an analogous coefficient matching for the tensors in Equation (60). The left-hand side of 
(60) is given by 

Q R ab{c d e) + Q K ab(c d e) = (X>i - V 3 + f s )A {a ri b)(c 5 e) d + (2? 3 -Th- f s )A^S b ^ d if e + - V 3 + f s )6^ dV b ^ c A^ 

+ 2f' u if b i 1 ce A d - 2/^ e( y )c I d - f sV ab A( c 6 e \ + (v 4 - J(x> 5 + v 7 )) A a A b A c s e) d - 2V 4 A { - a 5 b \A c A e 



l -V 5 - 2? 6 - f')j A a A b if £ A d - CD 5 - 2f' s )A^ c A^A d ~ Qz> 7 + A v ab A c A e A d - V 8 A a A b A c A e A d , (B5) 

and the right-hand side is given by 

(V 2 C 3 - 2V 3 )A^ a v b ^ c 8 e \ + (V2C1 - 2V 1 )A {a 5 b) d jf e + Vi(d + A 2 C 6 ) V ab v ce A d 

\v 2 C 5 - (X> 5 + V r ) + V 3 \C 3 + \A 2 C^j + U 3 \V 3 + ^A 2 (V 5 + 2? r )J J AM 6 I( c <5 e ) d 

- F x ^C 3 + \A 2 C^j ij^A^S^d + (V 2 C 4 - 2V 4 )A( a 5 b ) d A c A e + (^V 2 C 5 - V 5 - V^j A {a -q b){ - c A e) A d 

(v 2 c 6 - 2V 6 + v 3 (d + A 2 c 6 ) + u 3 {v t + A 2 v 6 ))A a ASf e A d + v x (c 4 + ^c 5 + A 2 c^j rf b A c A e A d 

v 2 c 8 - 2V 8 + v 3 (c 4 + ic 5 + A 2 c^j + u 3 (v 4 + l(v 5 + v 7 ) + A 2 v^j ^ A a A b A c A e A d . (B6) 

We have used the fact that both candidate vector kinetic terms found in the previous section have U\ — 0, U 2 = — 2 
and C 7 = 0. | 
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